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I. INTRODUCTION 



The neutrino masses and mixing which are required for giving a consistent explanation for 
the solar and atmospheric neutrino anomalies are the most firm evidence of physics beyond 
the electroweak standard model (ESM). New physics can be implemented in a variety of 
different scenarios. There are basically two main schemes that are often followed: (i) new 
matter content is added to the model respecting the original ESM gauge symmetry and (ii) to 
consider a model with a larger gauge symmetry. Certainly both schemes can be implemented 
together. In this vein, extensions of the ESM having an extra U{1) gauge symmetry factor 
are interesting for a variety of reasons. They are the simplest way of extending the ESM 
gauge symmetry and can be thought of as an intermediate energy scale symmetry coming 
from the breaking, at a higher energy scale, of a larger gauge symmetry describing some yet 
unknown physics. For instance, U{1) gauge factors are contained in grand unified theories, 
supersymmetric models, and left-right models. One major feature of these models is the 
existence of an extra neutral vector boson, usually denoted by Z', whose mass is related to 
the energy scale of the extra U{1) symmetry spontaneously broken. It is expected to have 
Z' signals at the TeV scale and its discovery is one of the goals of the LHC and future lepton 
colliders. Depending on the implementation of this kind of model, it can have a natural 
candidate for dark matter (DM) and/or furnish a mechanism for leptogenesis. Through the 
years much work has been done considering the features of this extra U(l) gauge factor and 

nn 

some particular formulations of the model were made. See, for example, Refs. |2|. In 
particular, when the charge of the extra U{1) factor is identified with B — L (baryon number 
minus lepton number), there is extensive literature concerning the most different versions of 
the model and a large variety of phenomenological aspects. 

In this paper we consider a B — L gauge model which has the particularity of being 
rendered anomaly free by introducing right-handed neutrinos with exotic B — L charges. The 
number of right-handed neutrinos and their B — L exotic charges is fixed by the anomaly 
cancellation equations. Since in this model not all of these right-handed neutrinos have the 
same exotic charge, we can construct Yukawa terms with different SU{2)l scalar doublets. 
Appropriate SU (2) l scalar singlets are also introduced to write the most general mass terms 
for the right-handed neutrinos. We make a detailed study of the scalar potential, concerning 
the mass spectra and the physical Goldstone bosons, and take advantage of this rich scalar 
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sector to construct a seesaw mechanism at low energies (TeV scale) to give realistic masses 
to the light active neutrinos. 

The outline of this paper is as follows. In Sec. II we present the particular B — L gauge 
model under consideration. In Sec. Ill we analyze the scalar potential of the model -the 
symmetries, the mass spectra and the model compatibility with experimental constraints 
-and introduce a Z2 symmetry to allow the model to have DM candidates. In Sec. IV 
we study the neutrino mass generation and show the compatibility of the model with the 
observed neutrino masses and the tribimaximal mixing. Finally, our conclusions are given 
in Sec. V. 

II. THE MODEL 

We consider the model of Ref. j3i] that we briefly summarize here. The model is an 
extension of the ESM based on the gauge symmetry SU{2)l ® U{1)y' ® U{1)b~l where B 
and L are the usual baryonic and leptonic numbers, respectively, and Y' is a new charge. The 
values of Y' are chosen to obtain the ESM hypercharge Y through the relation Y = [Y' + {B — 
L)], after the first spontaneous symmetry breaking. In order to make the model anomaly 
free we have to introduce right-handed neutrinos (n/j). Solving the anomaly equations we 
find that the number of riji cannot be less than 3, if we restrict ourselves to integer quantum 
numbers only. For the minimal number (3) these equations have two solutions: the usual 
one, where all right-handed neutrinos are identical and have L = 1, and the exotic one, where 
two of them have L = 4 and the third one has L = —5. The model under consideration has 
the right-handed neutrinos having such exotic lepton numbers. 

The fermionic content of the model is the same as the ESM plus the right-handed neu- 
trinos introduced above. The respective charge assignment is shown in Table HI In the 
framework of a gauge theory with spontaneous symmetry breaking, we at least have to 
introduce a scalar doublet, H, in order to give mass to the lighter massive neutral vector 
boson (Z) and the charged fermions, as in the ESM. However, more scalar fields are needed 
to give mass to the extra neutral vector boson {Z'), which is expected to be heavier than 
Z, and to the neutrinos of the model. Respecting gauge invariance, a general choice is to 
introduce two SU{2) scalar doublets, $i,2; and three SU{2) neutral scalar singlets, 01,2,3, 
with the charge assignments shown in Table HIl With these fields, and omitting summation 
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Table I: Quantum number assignment for the fermionic fields. 
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symbols, the most general Yukawa Lagrangian respecting the gauge invariance is given by 



£y = Y^'LLieRiH + Yl-'Q^idR^H + Y^Y'Ql.urjH + V,mLLinRm<l>i + VaLuUR^^i 



.(d)- 



+Mmn{nR^YnRn(t)l + A^33('^ra)'''^i?302 + A^m3('^i?m)'^'^ra03 + H.C. 



(1) 



where z, j = 1, 2, 3 are lepton family numbers and represent e, yU and r, respectively, m, n 
1, 2, and if = iT2H* . The corresponding scalar potential is 



Vb- 



B-L 



-filH^H + Xu^H^Hf - + An 



/i22$2'^2 + A22 



01. I V^Q/'rQ: I 



A12 |$l|' |$2|' + A;2($I$2)($5'^'i) + ^H, \H\^ 



^ + A' 1$ P 



+ /3l230102(03)' + H.C.] + A,^(0>,)(0;0;3) 



+ 



$I$2(/3l3010; + /3230;03) 



(2) 



where 7 = 1,2; = 1,2,3; and a < /3 in the last term. In Cyi the motivation for 

introducing such scalar fields is to write the most general neutrino mass terms. Because of the 
fact that not all right-handed neutrinos have the same Y' and {B — L) charges, the neutrino 
mass matrix will have entries proportional to vacuum expectation values (VEVs) which can, 
in principle, belong to different energy scales. The scalar potential is a consequence of the 
fields we have previously introduced, and the terms in Eq. ([2]) are only dictated by gauge 
invariance. Now, we have to observe that when we write terms based on general grounds, 



although correct, we may have introduced more symmetries than we need. Hence, we must 
do a detailed study of the scalar potential, and know the scalar mass spectra in order to 
avoid inconsistencies with the present phenomenology. The scalar doublets of the model, H 



Table II: Quantum number assignment for the scalar fields. 

h I Q Y' B-L Y 

=Fl/2 1/2 0,1 1 1 

^>°" ±1/2 1/2 0,-1 -4 +3 -1 

$2' ±1/2 1/2 0,-1 5 -6 -1 

01 -8 +8 

02 10 -10 



0001 -1 



and $1,2; contribute to the Z boson mass, so their vacuum expectation values are bounded 
by the electroweak energy scale. Hence, the largest energy scale of the model comes from the 
SU{2) scalar singlets 0i,2,3- In this way, the pattern of the spontaneous symmetry breaking 
is 

SU{2)L(i^U{l)y,®U{l)B-L (01.2.3) ^ 

SU{2)l^U{1)y (i/,$1.2) f/(l)em. (3) 



III. THE SCALAR POTENTIAL ANALYSIS 

Now, we focus on the analysis of the Vb-l scalar potential given in Eq. ([2]) when all neutral 
scalar fields develop nonvanishing VEVs, with the usual shifting = -^{V^ + Keip + ilmip). 
By using standard procedures we are able to find the constraint equations coming from the 
linear terms in the scalar potential after the symmetry breaking. See the appendix. In the 
same way, we can construct the mass-squared matrices for the charged, real, and imaginary 
scalar fields. We start looking at the mass-squared matrix for the charged fields. It is a 
complete 3x3 symmetric matrix in the basis {H^, <l>^) that can be easily diagonalized 
and, after taking into account the constraint equations, gives the following mass spectrum: 
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two charged Goldstone bosons 

y *2 *2 

and two massive states whose expressions we are not showing by shortness. The fields 
will be absorbed to form the longitudinal components of the charged massive vector bosons 
W^. The other two physical states remain in the spectrum and are a prediction of the 
model. Later in the paper we approach numerically all the mass spectra in some different 
situations. 

In the neutral imaginary scalar sector we have a 6 x 6 mass-squared matrix that, after 
the diagonalization procedure, shows two massive scalar and four massless fields. Two of 
them will become the longitudinal components of the Z and the Z' neutral vector bosons. 
The other two massless states remain in the physical spectrum. We show the two physical 
Goldstone bosons only in the limit where V^i ^.s ^ Vh, V^i,2 they are given by 

G'^ ^ . / t . ( + 4^^Im$? + 4^^Im<|.°l , (5) 



G^p^ ^ ^=(7Im0i + 5Im02 + 6Im03). (6) 
V 110 

From the expressions above we see that G%^ ^ are mainly doublet and singlet, respectively. 
This fact will be analyzed later on. 

For the neutral real scalar sector we have a symmetric 6x6 mass-squared matrix with 
a nonvanishing determinant. Hence the spectrum will not contain massless states. We will 
not show analytical expressions here but we will give numerical values below. 

The number of Goldstone bosons we have found by doing explicit calculations can be eas- 
ily understood by studying the global symmetries of the scalar potential before and after the 
spontaneous symmetry breaking (SSB). Before the SSB, the global symmetries of the scalar 
potential are (a) SU{2) acting on H and $12 doublets, (b) U{1) acting on H with charge +1, 
(c) U{1) acting on <l>i 2 with charge +1, and (d) two independent U{l)j3^^ transformations 
acting on the fields $1, $2, 0i, 02, 03 with charges (|, — |, 1, —1, 0) and (— |, |, 0, +2, 1), re- 
spectively. After the SSB the global symmetries of the scalar potential are reduced to a single 
U {l)a acting on the charged components of the doublets, H"^ and $^2 with charges (±1, ±1). 
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Following the Goldstone theorem, the number of Goldstone bosons is equal to the number 
of broken symmetry generators. In this case the original symmetry has (3 + 4 x 1) = 7 gen- 
erators and the remaining symmetry has 1. Then we must have 6 Goldstone bosons, which 
is exactly the number we have found just above: two charged and four neutral imaginary 
fields. 

Notice that the scalar potential given in Eq. ([2]) corrects the one given in Eq. (16) of 
Ref. in which the terms proportional to A'^:^^ are missing. The lack of these terms alters 
the global symmetries under which the scalar potential is invariant and, consequently, the 
number of Goldstone bosons in the spectra. In that case the symmetries before the SSB 
are (a) 0(4) acting on the four components of (b) SU(2) acting on $i,2) (c) U{1) acting 
on $1,2; (d) the two [/(l)^ ,^, defined above. After the SSB the remaining symmetries are 
(i) 0(3) acting on the components (Imif'^, KeH^, ImH^) and (ii) U{1) acting on $^''2 with 
charge ±1. Therefore, we are left with (6 + 3 + 3x 1) — (3 + 1) = 12 — 4 = 8 Goldstone 
bosons. The same result is obtained by doing explicit calculations. From the mass-squared 
matrices we find that the number of Goldstone bosons is the expected one and also that in 
the charged sector we are left with four massless states given by 

= H^, (7) 
= ^ {V^M + V^M) ■ (8) 

^yi + yi 

Hence, this result is in conflict with the present phenomenology since there are two extra 
charged massless scalars in the spectrum. 

Now, we return to our present analysis. Since our analysis of the scalar potential shows 
the existence of two physical Goldstone bosons, it is time to care about the safety of the 
model. Before that, some remarks about the VEVs of the model are due. The V^^ ^ 3 are the 
largest energy scale of the model. The main contribution to the Z boson square mass comes 
from the doublets so that (V^ + Vi^ + V^J = Vi^u = (246)^ GeV^ The doublet H is the 
one that couples to quarks and to charged leptons via Yukawa interactions, and hence, Vh 
must be close to Vesm to give the correct tree level mass to the quark top, as the ESM do, 
for an 0{1) top Yukawa coupling. We then conclude that V^^ + V^^ <^ V^. 

The major challenge to models with physical Goldstone bosons, also called Majorons (J), 
comes from the energy loss in stars through the processes 7 + e~ — )■ e~ + J. This process is 
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used to put limits on the eeJ coupling, and it is found that it has to be g^ej < 10"^" for the 
Sun, and g^ej < 10^^^ for the red-giant stars [4]. However, the dynamics of the red giants 
has not the same level of confidence as that of the Sun, and this fact considerably weakens 
the second constraint. 

The physical Goldstone has components only in the SU{2) singlets 0i,2,3, which 
couple only to right-handed neutrinos. Therefore, it is safe since there is no tree level 
contribution to the energy loss process. The case for G^p^ is not that simple. has 
a component in the ESM-like doublet H, and it contributes to eeJ through ImH^. The 
components in $1,2, which couple only to neutrinos at the tree level, pose no problem. 
Since in this case symmetry eigenstates and mass eigenstates are connected by orthogonal 
matrices, from Eq. ([6]) we find 



\/2 J- T/2 

ImH^^J , , "^' Gl + ..., (9) 



and hence, 



^ X^J K^i. . ( ^ ] ^ 2 X 10-^ < 10- - 10-". (10) 



^\vli + vl + vl^ \ V„ l\V2V„J v„ 

where is the electron Yukawa coupling to the doublet H, we have defined Vj^ -|- V^^ = V^, 
and we have used the shift H'^ — )■ + Reif° -|- ilmH^). From the equation above we 

conclude that the VEVs of the SU (2) doublets $1^2 must be less than 12 MeV to satisfy the 
Sun constraint or less than 120 KeV to satisfy rigorously the red-giant constraint. Let us 
adopt for practical purposes an intermediate scale: V$ = 1 MeV. 

Once we have established the energy scale of the VEVs of the model, and verified its 
safety up to now, we can make an exemplary study of the full scalar mass spectra. We will 
do it numerically since the excessive length of the analytical expressions make them useless. 

In order to compute the masses we consider a set of parameters: the dimensional ones 
Vh = 246, V^,^, = 0.001, V<^,,,,3 = 1000, in GeV, and = 0.2, An,22 = A'i3,2i,22,hi,//2 = ^sa = 
1, ^Hi,H2 = ^Hsa = A'^^ i2,23' ^12,13,23 = 0.1, Pus = "0.8, which are dimensionless, for 
a = 1,2,3. Note that the values for the fi parameters are found by solving the constraint 
equations for the scalar potential: fi'j^ = (402)^, fif-^ = —(630)^, /i22 = (230)^, fi^i = /i^2 = 
(838)2, ^2^ = (550)2, in Q^y2 ^^^^ ^ ^ 0.4885, and g = 0.6298, where gy, 

dB-Li and g are the coupling constants of the t/(l)y', U{1)b-Li and SU(2)l gauge factors, 
respectively, which are related to the electric charge through 1/6^ = l/gY' + i/gB-L~^^/9^ 
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The charged scalar sector gives the masses mc^. = (1424.9, 173.9, 0), also in GeV, where 
the massless complex field is responsible for the longitudinal components of the charged 
vector bosons W'^ and W' . 

In the imaginary neutral scalar sector we have a 6 x 6 square mass matrix and after 
diagonalization we find, in GeV, rrij^ = (1549.2, 1414.13, 0, 0, 0, 0). Notice the correct 
number of the massless fields: two are absorbed to form the longitudinal component of the 
neutral vector bosons of the model, Z and Z', and the other two are the physical Goldstone 
bosons G^p^^, as discussed above. 

In the real neutral scalar sector, in the same way, we find m^j. = 
(1743.9, 1643.2, 1414.2, 1029.8, 150.0, 0.0014), in GeV. We have found a very light 
scalar of about 1.4 MeV which poses a new challenge to the model: the Z invisible decay 
width. The presence of such a light scalar field, say R, and the G%_^ = J zero mass state, 
allows the decay Z RJ — )■ J J J, which will contribute to the Z invisible decay width as 
half of the decay Z — )■ uu, for a single fiavor family j^l. According to the experimental data 
there is no room for such an extra contribution p]. 

The light scalar we found above is not the result of a particular choice of the input param- 
eters, as it could be thought at the first moment. Let us provide a qualitative but convincing 
argument. As was observed in Ref. ^], the reason is as follows. We have mentioned above 
that, before the SSB, the scalar potential has a U{1) global symmetry acting on each of the 
$1.2 doublets, say, $. This means that we can rotate freely in the Re^^-Im^*^ plane, so 
that as long as this f/(l) symmetry holds, the fields Re$° and Im^'' are mass degenerate. 
However, this symmetry is broken when the real neutral component acquires a nonvanishing 
VEV and, hence, the fields are no longer mass degenerate. The square mass difference must 
be, then, of the order of the square of the energy scale responsible for breaking the symme- 
try, i.e., ■m^g^o — ""^im.j.o = ^(yi)- When Im$° is a Goldstone, mf^^o = 0, we are left with 
m^g^o = C>{V^), which, in our case, it is a very light scalar since V<j, must be of the order of 
1 MeV, in order to be consistent with the star energy loss data. Then, we must find a way 
to reconcile the present model with the experimental constraints. 

Some attempts can be made to remove such inconsistency. Since the origin of the problem 
is in the breaking of the U{1) symmetry acting on the doublets $1,2, let us consider the 
situation where V$2 = 0, and all other VEVs are different from zero. In this case we find 
the same number of neutral Goldstone bosons (4): two would be Goldstone bosons and two 
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physical ones G'p_^ ^ . G'p^ is given by the same expression as in Eq. ([6]) , and 

G"^ ^ , ^ (I4,lmiy° + Vnlm^T) ■ (H) 

We also find that for the same input parameters, but now providing an input value for 
/^22 ~ (230)2 GeV^, the mass spectra are practically not affected and we still have a light 
real scalar whose mass is about 1 MeV~ (9(V$J. We get the same conclusion if we consider 
= and V$2 ^ 0. We only have to do the replacement $i < — > $2 in the above results. 

As the problem persists, let us now consider the case where = and V$2 = 0. 
In this case, the number of Goldstone bosons is reduced to 3. There is only one phys- 
ical Goldstone, the given in Eq. (Q, which is safe, as discussed above. The mass 
spectra are now considerably affected. For the same input parameters as above, and 
with fi^-^ = —(800)2,^^22 = (230)2, in GeV^, the spectra, with all the masses in GeV, 
are the following. For the charged scalars we have rric^ = (1469.4, 380.1, 0). For the 
imaginary scalars we find rrij^ = (1549.2, 1459.1, 337.9, 0, 0, 0), and for the real scalars 
rriR^ = (1743.9, 1643.2, 1459.1, 1029.8, 337.8, 150.0). As we can see, there is no a light real 
scalar anymore. The lighter real scalar is heavier than the Z vector boson, so that the prob- 
lematic decay Z — ?► RJ is kinetically forbidden. Then, we have succeed in making the model 
safe. However, this solution is not satisfactory since the choice we have made for the doublet 
VEVs (V$j 2 = 0) does not allow the light neutrinos to get mass. It is easy to see that in 
this case there is a remaining U{1) quantum symmetry, say, f/(l)f, protecting the neutrino 
mass generation at any level. A possible (^-charge assignment is: C(^^eL; e^, C/j, $1^2) = ~1, 
C{uL,dL,UR,dR) = 1/3, and C(^(i,2,3)r, 0i,2,3) = 0. In order to make the model compatible 
with the experimental data and, hence, with massive neutrinos, we have to look for a new 
kind of solution since the symmetry breaking pattern above is not realistic. 

Before continuing the search for a satisfactory solution, let us observe that before the SSB 
the model has a Z2 exact symmetry with the transformation rules ^2(71/23) = —riji^, Z2($2) = 
— $2) ^2(^3) = ~03; and all the other fields being even under Z2. It is interesting to 
preserve this symmetry after the SSB if we are looking for DM candidates. This is true 
when V$2 = V^a =0. In this case the Z2 symmetry is not spontaneously broken, and a 
mechanism similar to that of Ref. jll can be implemented. The number of Goldstone bosons 
is 4, and the physical ones are given by the following: is given by the same expression 
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in Eq. ([TT]), and 

G% = . ^ {W^M<Pi + 4l^^,Im02) . (12) 

However, as we already know, there is a very light real scalar that, together with G%-^, has 
severe implications on the Z invisible decay width. We will come to this Z2 picture later on. 

A. The solution 

With the aim of constructing a consistent model, let us introduce a new SU (2) neutral 
scalar singlet (f)x with the quantum numbers Y' = —{B — L) = ?>. The Yukawa Lagrangian 
remains as in Eq. ([T]), but to the scalar potential in Eq. ([2]), besides extending the range of 
the indices to a, /9 = 1, 2, 3, X, we have to add the following non-Hermitian terms, 

V^_L = -^^^Hlx^lr2H<Px - iKH2x{^lr2H){<P\f + /3x(0x0i)(0203) + P^x{<P*x<Pl) + H.c, 

(13) 

in order to account for all the gauge invariant terms after the introduction of (px- The terms 
above reduce the number of global symmetries of the scalar potential, so that changes in 
the scalar spectra are expected. 

Before the SSB, the global symmetries of the total scalar potential are (a) SU{2) acting 
on H and $12 doublets, (b) U{l)a acting on if and $1^2; and (c) f/(l)/3 acting on the fields 
if, $1, $2, 01, 02, 03, 0x with charges (|, 0, — |, 1, — |, — |, — |), respectively. After the SSB 
the global symmetries of the scalar potential are reduced to a single U{1) acting on the 
charged components of the doublets, and $^25 with charges (±1, ±1). The total number 
of Goldstone bosons will be given by the number of broken generators, i.e., 5 — 1 = 4, which 
is the number of massless fields needed to form the longitudinal components of the charged 
{W^, W~) and neutral vector bosons {Z, Z'). In this case, there are no physical Goldstone 
bosons at all. It means that the inclusion of the SU{2) scalar singlet (px has removed all 
physical massless states from the spectrum, and we have succeeded in finding a solution for 
the safety of the model. 

Now, numerical applications require expanding the input parameters set to account for 
the new ones related to (px- We then choose V"^^ = 1000, andnHix = 0.01 in GeV, 
and the dimensionless Xsx = A'lx = 1, ^HsX = Kx = ^ix = ^2X = ^3X = f^sx = 
0.1, f3x = — 0.6, and ftiy2X = 0.001. As before the /x parameters are found by solving the 
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constraint equations given in the appendix. With the above parameter set, plus the one 
we have used previously, we find mc^ = (11 137.3, 1661.7, 0) for the charged scalar sector, 
rrij- = (11135.9, 1652.6, 1467.0, 973.6, 0.002, 0, 0) for the neutral imaginary sector, and 
niR^ = (11 135.9, 1927.6, 1816.6, 1652.7, 1508.8, 900.5, 146.2) for the real scalar sector, in 
GeV. Notice that we have now a very light pseudoscalar, which has components mainly in 
the SU{2) singlet fields 0i,2,3,x- For instance, its component in ImH^ is 7.3 x 10^^^, which 
implies Qeej ~ 10~^®. Hence, it is compatible with the astrophysical constraint, and poses 
no problem to the Z invisible decay width, since all the real scalar fields are heavier than 
the Z boson. 

In this case the introduction of the (px scalar provides the right elements to make the 
model safe. Moreover, concerning the neutrino mass generation, from the Yukawa terms in 
Eq. ([1]) we are able to construct the most general neutrino mass matrix, since now all VEVs 
are different from zero. 

B. A Z2 symmetry and dark matter 

Now let us consider the Z2 symmetry again, after the introduction of the scalar (px- We 
had the field symmetry transformation rules 

Z2{nR3) = -nm, ^2(^2) = -$2, ^2(03) = -03, (14) 

and all the other being even. It is easy to see that all the Hermitian terms in the scalar 
potential involving (px are invariant under Z2. However, the non-Hermitian terms 

-ZKH2x{^^r2H){4>*^Y, /3x(0x0i)(0203), and f3M(px<Pl), (15) 

in V^_i^, are not invariant. We could change the (px transformation rule to odd, in order 
to have some of them invariant. In this case, however, if we want to keep the Lagrangian 
invariant under Z2 after the SSB, we must have V<^^ = 0, and this is not an option since we 
need V^^ 7^ to have a consistent model, as discussed above. Motivated by the possibility 
of having DM candidates we impose the Z2 symmetry to the entire Lagrangian. Then, the 
terms in Eq. (fTSj) will be removed from the scalar potential and the only non-Hermitian 
terms allowed are 

y^^.i = $l$2(/3l30103 + Mlh) + ^230102(0;)' - tKHlX^^T2H<Px + H.C. (16) 
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After the SSB, the Z2 symmetry is not broken if we have V$2 = V^g = 0, and we have mass 
eigenstates that are also eigenstates of this symmetry. However, we know from our previous 
analysis, before introducing 0x, that this vacuum configuration challenges the safety of the 
model due to a physical Goldstone and a light real scalar. Now, after introducing 0x we 
have four massless states in the neutral imaginary sector. However, in this case, both of the 
physical massless states are mainly singlets: is given by the same expression in Eq. (fT2l) . 
and 

G% ^ Jl (12 Im0i - 15 Im02 + 82 Im0x) • (17) 

' V7093 

In fact, for the parameter set we used above, the G%_^ component in ImH^ is ^ 2x 10^^^ which 
implies Qeej ^ 4 x 10^^^; thus it is safe with respect to the star energy loss constraint. This 
main feature is due to the introduction of the trilinear term —iK,Hix^i'T'2H(j)x- Qualitative 
arguments to explain this behavior can be given. The number of U{1) symmetries is the 
same in both situations, with and without (px, 4. Without (px, we have two independent 
f/(l) symmetries involving only the doublets, say, f/(l)o- acting on H, and U{l)a acting on 
$1.2. With (px, the trilinear term relates the U{1) charge of H to that of $1, reducing the 
number of U{1) symmetries involving only doublets to just one, say, U{l)a acting on H 
and $1,2, and at the same time, it introduces a new U{1) symmetry acting on H and (px, 
say, U{l).y. The number of broken generators is the same in both situations, since we have 
the same number of massless states; however, the origin of these physical massless states is 
different. The introduction of (px works in a very similar way to the singlet introduced in 
Refs. [8] and [9] to form the terms H'^T2Hd(p and H^T2Hd(p'^, respectively, in order to make 
the axion invisible. 

The numerical spectra for the different scalar sectors, in GeV, 
are mc^ = (1489.9, 1330.3, 0) for the charged scalar sector; m/,. = 
(1479.7, 1433.9, 1318.9, 0, 0, 0, 0) for the neutral imaginary sector, and m/j. = 
(1483.5, 1479.7, 1378.4, 1378.4, 1318.9, 675.3, 152.9) for the real scalar sector. The 
point here is that all real scalar fields are now heavier than the Z boson, avoiding in this 
way the Z invisible decay width constraint. Therefore, the model is safe from these most 
severe constraints. 

As the Z2 symmetry still holds after the SSB, due to this particular vacuum configuration, 
the model can present some DM candidates. In general, a candidate must be the lightest 
particle odd under Z2, in order to be stable. In our case, it can be the lightest odd mass 

13 



eigenstate of the nns or the hghtest odd imaginary mass eigenstate, or its odd real counter- 
part. This subject will be considered elsewhere. Here we only estimate the relic abundance 
and the direct detection cross section for the case of the fermionic cold DM candidate ur^ 
(referred as x from now on). 

The most relevant annihilation process of x occurs via the t-channel exchange of $^($2) to 
charged (neutral) leptons' final states. The thermally averaged x annihilation cross section, 
(cTf), is given by 10| 

{av) ^a + b (v^) ^^Y: ^e.,.- K ■ (18) 

■ij 

where i,j = e, /i, r and Cc are the color factors, equal to 1 for leptons. Also we have neglected 
the lepton masses. The GefT.y = Di^D*.^/ {m^^ + ^x) effective couplings, where 

mc^ ~ 1489 GeV {mn^ 1489 GeV), the mass to be considered when c harg ed (neutral) 



leptons are produced. The relic abundance of x is approximately given by 

^,2 1.04xl0-9xf 



(19) 



Mpiy/g^{a + Sb/xf) ' 

where, in this model, g^, = 107.75 is the number of relativistic degrees of freedom available 
at the freeze-out temperature, Ty, and Xj = M^/Tj is given by 

fA5g^MpiM^{a + 6b/xf) 



In 



(20) 



with c = 5/4 and = 2. Using the following set of parameters, D^s = 0.06, = 
0.9, De3 = 1, and for = 750 GeV, we find xj = 24.81 and ilh'^ = 0.11, which is in agree 



ment with the experimental bounds 



121 . The same interaction allowing the x annihilation 



in charged leptons, which are proportional Di^, also induces lepton flavor violation (LFV) 
such as yU — )■ 67 and t ^ fi'y (see below). 

The next task is to compute the direct detection cross section. In our case, the elastic 
scattering of x with nuclei occurs via the t-channel x + — )■ x + process due to the 
exchange of the scalar mass eigenstate Rj, which is the Higgs scalar boson in the model with 
mass m/j^ ^ 152.9 GeV. The spin-independent cross section is given by |13| 



XN 



[Zfp + {A-Z)U 



vr (M^ + m^y 

where the effective couplings to protons and neutrons, are 

J To ~ r '^JTG 2-^ 



fp,i 



E 

q=u,d,s 



q=c,b,t 



a/2 ruq 



(21) 



(22) 
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In this case Gefr.g = Go x = [C^j^iRiCHR^Mj (V^j V^m^J] x m^, where C^^r^ ^ 0.01 and 
Chrj ~ 0.99 are the coefficients relating the symmetry eigenstates (02? H) to the relevant 
mass eigenstate i?7, respectively. By using f!^^"* and Z^^^"'* given in Ref. y^] we find 



3 X 10"^ pb X 



Geff,g X (1 GeV/mq 



n 2 



(23) 



which gives cTy^p ^ 4.74 x 



recent present bounds 



10-7GeV-2 

0~^^ pb, for = 750 GeV, which is in agreement with the most 
15Nl7|. The parameter set we have used in all the cases above is 
compatible with the following requirements: (i) the constraint equations are satisfied, (ii) 
all obtained masses are iv? > 0, and (iii) results for the already known fields are consistent 
with those of the SM at the tree level. 

Notice that, although we are considering a multi-Higgs model, the values we have found 
for the lightest real scalar, the Higgs boson, are in agreement with the last combined CDF 
and DO results for the ESM Higgs boson, which have excluded, at the 95% C.L., a region at 



high mass in 158 < m^/ < 175 GeV 
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IV. NEUTRINO MASSES 



The model without the Z2 symmetry already has a satisfactory solution to the neutrino 
masses, since we are able to construct a general neutrino mass matrix. However, we are 
going to consider the case with this symmetry because the model becomes more attractive 
due to the presence of stable candidates to DM. 

The Yukawa Lagrangian in Eq. ([T]) gives the following neutrino mass terms: 



- Cm, = Viml^LinRmV^^ + Mmn{n'^m) LriRnV^^ + Mzz{nl) lUr^V^^ + H.C., (24) 

where z,j = 1,2,3 (or e, /x, r, respectively, when convenient) and m,n = 1,2. In matrix 
form Eq. reads 



_ r 

'r; 



Mr) 

M 



Ml M, 



riR 



with 



ur = [rii n2 ns]^ . 



(25) 



(26) 
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The Majorana mass matrix (Mm) and the Dirac mass matrix (Mr)) are given by 

^ Mn Mi2 ^ 

Mm = V^, Mi2 M22 , Md = V^ 



\ 











M 



33 



^ Vn Vi2 ^ 



#1 



V21 V22 

I?31 I?32 



(27) 



since V^^ = = 0, where Mm = Mj^. For ^ V^/.^, the mass matrix in Eq. can be 
diagonahzed by an approximate scheme. The masses of the heavy neutrinos are related to 
the energy scale of the VEVs of the singlets 0i and 02? and are given by the eigenvalues of 
Mm: 

{Mu + M22) T V4A^?2 + (-Mil --^22)2. 



Ml 



M. = M,,,v: 



The masses of the light neutrinos are given by the eigenvalues of the matrix 
which are 



1 



2D 



M 



A T VA2 + 



where the following definitions have been used: 
^1 = (Pn, P21, ^^3i), "^2 = (^^12, ^^22, ^^32), r = ADm 

A = Mii(^2f + A^220i)' - 2Mi2(^1^2)) D, 
Dd = ^^21^^32 — ^^22^^31- 



ms = 0, 



(^^12^1 



V 



11 



-M?2 - -^11-^22, 



(28) 



(29) 



(30) 



The parameters in Mm and M^, have to be chosen in order to have light neutrino masses 
consistent with the solar and atmospheric experimental data. However, since there is no a 
standard procedure to do that, we will present a particular solution to show that this model 
can generate realistic active neutrino masses. 



A particular solution 



From the experimental neutrino data it is found ^ that the neutrino mixing matrix is 
compatible with the so-called tribimaximal (TB) one [19j, which is given by 



TB 



2 J_ 

3 







\ 



^_ j_ ]_ 

Ve Vs V2 

. L J_ A. , 



(31) 
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and where it is assumed that the neutrino mixing angles are in a good approximation given 
by sin^^ia = 1/3, sin^ ^23 = 1/2, and sin^ ^13 = 0. Working in a basis where the charged 
lepton mass matrix is aheady diagonal, the U^b matrix diagonalizes the light neutrino mass 
matrix in Eq. (125]) : U^^M,^Utb = My = diag(mi, m2, ms). It can be shown that the most 
general neutrino mass matrix that can be exactly diagonalized by ?7tb has the form 



Mr 



TB 



^ X y y 

y X + u y — V 
yy y — I' x + u 



\ 



using the same notation as in Ref. 
The Mtb mass eigenstates are 

mi = X — y, 



19| 



m2 = X + 2y, 



= 2u + X — y. 



(32) 



(33) 



The square mass differences Am^Qi and Am^^-j^^, needed to explain de solar and atmo- 
spheric neutrino anomalies, can be obtained by imposing conditions on x, and v. The 
simplest way to apply this analysis to our particular case is as follows. We consider 

V^2 



Mu = M 



22 



-M 



33, 



M 



12 



0, 



so that the Majorana and Dirac mass matrices are now given by 



Mm = MuV^,1 



3x3j 



D 



( ^ 

V21 V22 
V^i V32 



(34) 



(35) 



Then, the light neutrino mass matrix becomes 



\/2 1 



V^, Mn 



VI 



VnV2l + Vi2V22 Vi, + ©22 



VnV2i + VuV22 VnV^i + VUV32 



VnV^i + V12V32 V21V31 + 1^22^^, 



32 



J 



(36) 



The matrix above has a null determinant and, therefore, a zero mass eigenstate. Hence, in 
order to make both matrices compatible, we must have a vanishing eigenvalue in Eq. (133|) . 
We choose my, = 0, i.e., 2^ + x — y = 0, and, hence, x + v = y — v. 
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Comparing Eq. ( !32|) with Eq. (136|) we have the following equations: 

^ = ^11+ ^'u, (37) 
^ VnV2i + V12V22 = VuV-ii + 1^122^32, (38) 



K 



-Vl + Vl, = Vl,+Vl,, (39) 



K 



V^lV^l + 2^222^32, (40) 



K 



where we have defined the dimensional constant K = . . ,V . A solution for the above 



equations is 

V2i = Vsi, and V22 = Vs2. (41) 

From the above equations we see that the condition to have = 0, x + u = y — u, is 
automatically satisfied. We have the following equations to fit the atmospheric and solar 
neutrino data, 

mi = X — y, 1712 = X + 2?/, m-^ = 0, (42) 

and, therefore. 

Ami, = ml-m\ = 3y{2x + y) > 0, (43) 

l^^ltml = = -y)'^- (44) 

Assuming that x — y > we have to solve the equations 

3y{2x + y) = 7.67 x 10"^ (eV)^ and x-y = (2.4 x 10"^)^/^ eV, (45) 

which are satisfied by x = 0.049 248 7 and y = 0.000 258 887, in eV. The corresponding 
mass eigenvalues are then given by mi = 0.048 989 8, m2 = 0.049 766 5, and m^ = 0, in 
eV, showing an inverse hierarchy pattern. We can now solve Eqs. (137|) and (138|) for the 
Vij parameters. In order to do that we have to know the value of the dimensional con- 
stant K. For = 1 MeV, V^^ = 1 TeV, and assuming Ain = 1, we have K = 1 eV. 
Choosing the input values V22 = 0.25 and V21 = 0.15, we find Vu = 0.190 751, and 
V12 = —0.113 415. Experiments on Oz//3/3 can put bounds on |mee|, and the strongest one is 
\mee\ < 0.26 (0.34) eV at 68% (90%) C.L. [20]. This quantity is related to the mass eigen- 
values through \mee\ = |cf 3(777.10^2^*'^^ + 777.25^2^**^^) + TTisc^*'^'^^ s^g | . In our case, with no 
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CP violation nor phases in the leptonic mixing matrix, we find |mee| ~ 0.05 eV. Future 
experiments, however, expect to improve sensitivity up to ^ 0.01 eV [21]. 

The procedure we have followed for finding a particular solution for the light neutrino 
masses can also be realized by using, instead of the matrices given in Eqs. (IHT]) and (152]) . the 



ones given in Ref. j22|, provided we make, in the notation of this reference, c = — (i/2, and 
the identifications u = d — {a + b) , y = d, x = a + 2b — d. It results —a = x — y + 2v = m^, 
and we take a=0. 

The results showed above demonstrate that the model is fully compatible with the exper- 
imental neutrino data, and that light neutrino masses can be generated neither appealing 
for very large energy scales nor imposing fine-tuning. Now, we have to verify if the set of 
parameters we have used above is in agreement with the LFV constraints coming from a 
process like li — )■ Ij + 7, where i = 2,3 = fi,T and j = 1, 2 = e, /i, respectively. This model 
has one loop contributions to such a process since charged leptons couple to charged scalars 
and right-handed heavy neutrinos. The branching ratio is estimated as j23| 

B {k ^ + 7) = (I/miI' + , (46) 

where a ~ 1/137 and Gp — 1-16 x 10~^ GeV~^ is the Fermi constant and 



t 4 (4^) 



with F2 (x) being 



^ , , I — <ox + ^x"^ + 2x^ — Qx'^lxix , . 

^^<^>^ ■ '''' 

Using the parameters needed to fit the neutrino masses and the ones to estimate Vth^ ^0.11 
(Pe3 ~ 0.06, V^,^ ~ 0.9, ^ 1, m„^3 = 750 GeV, mc^ ~ m^± = 1.33 TeV, mc^ = m^± = 
1.48 TeV), we can give an estimate for the branching ratio i? (/i — e -|- 7) ^ 7.9 x 10~^^ and 
-B (r — 7- /i -|- 7) ~ 2.5 X 10^^. These values are in agreement with the present upper bounds 



5 (/i ^ e + 7) < 1.2 X 10-11 and 5 (r ^ /i + 7) ~ 6.8 x 10^^ j^, 1251. 

The ratio between the VEVs we have used for finding the neutrino mass eigenvalues is 
V$^/V<^^ = 10"^. This is of the same order as the ratio mg/mtop = Y^/Ytop ~ 10^^, and it is 
comparable with rriu/mtop = Yu/Ytop ~ 10~^. We have chosen those values for and V^^ in 
order to have light neutrino masses without resorting to very tiny neutrino Yukawa coupling 
constants, or fine-tuning, and, at the same time, to have the Z' vector boson not extremely 
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heavy. This is a kind of seesaw mechanism where the heaviest scale, V^^, is constrained 
by the Z' vector boson, which should be not too heavy in order to not decouple from the 
spectrum. The light scale, , is then used to fix the absolute neutrino mass scale through 
the ratio V^^jV^^. In this picture we are substituting a hierarchy in the VEVs, which would 
have a possible explanation based on the dynamics of the fields, for one in the Yukawa 
coupling constants, for which we cannot find any natural explanation. This is basically the 
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philosophy behind the work in Refs. 

As we have discussed above, the absolute neutrino mass scale depends on the ratio V^^jV^^ 
where V^^ is a tiny value. Although this value can be affected by radiative corrections, it can 
be argued that, when the symmetry is considered, setting V$j to a tiny value, at the tree 
level, is natural because if it were in fact taken to be zero this would increase the symmetry 
of the entire Lagrangian ('t Hooft's principle of naturalness). This can be seen considering 
the constraint equations with V$j — > 0. It implies that Hu\x = 0, since Vh and V^^ 
differ from zero. Then the term —iKHix^iT2H(j)x does not appear in the scalar potential, 
Eq. ( fT6l) . and the entire Z2 invariant Lagrangian is now invariant under an additional global 
quantum symmetry, say, f/(l)(- A possible (^-charge assignment is (^(z/eL, e^, e/j, $1,2) = — 1, 
({uL,dL,Uji,dji) = 1/3, and C('^(i,2,3)R5 01,2,3) = 0. Thus, it is expected that the VEV 
hierarchy will remain stable when radiative corrections are taken into account. 



V. CONCLUSIONS 



In this paper we have studied in detail the scalar and the neutrino Yukawa sectors of 
an extension of the electroweak standard model which has an extra f/(l) gauge factor, as 
described in Sec. II. We have analyzed the scalar spectra of the potential given in Eq. ([2]) 
and found that it is inconsistent with the experimental data coming from the star energy 
loss and the Z invisible decay width. We would like to stress that this is a general result for 
this scalar potential. 

We find that the more suitable solution to this problem is the addition of a new SU(2) 
scalar singlet, called (px in the text. The new terms introduced by (px are able to remove all 
the physical Goldstone bosons and, at the same time, to have all the real mass eigenstates 
heavier than the Z boson. This solution is particularly interesting since, in this case, all 
VEVs can be different from zero, which allows for the construction of a general neutrino 
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mass matrix. 

In order to have a still more attractive model we consider the possibility of having DM 
candidates by including a symmetry. Before the SSB the only fields having odd trans- 
formation under are n/j3, $2? and 03. Z^ will still be a symmetry if the scalar fields 
$2 ai^d 03 do not develop VEVs. Hence, after the SSB we will have states which are mass 
and Z2 eigenstates simultaneously. It opens the possibility of having DM fields since the 
lighter Z2 odd eigenstate will be stable. Moreover, we show in a preliminary study that the 
fermionic field n^a is a viable cold DM candidate. 

We consider in detail the neutrino mass generation in the framework of the model with 
the Z2 symmetry. In this case we found an inverted hierarchy compatible with the solar and 
atmospheric neutrino data and the tribimaximal mixing matrix. Two appealing features 
are (i) the absolute scale of the neutrino masses is obtained by a seesaw mechanism at 
O (TeV) energy scale, which is the scale of the first symmetry breaking, and (ii) the observed 
mass-squared differences are obtained without resorting to fine-tuning the neutrino Yukawa 
couplings. 

The model has also some phenomenological implications. One of them is the existence 
of an extra neutral vector boson, Z\ which can be in principle detected at the LHC or 
International Linear Collider. In fact, there are studies showing that the Z' of this partic- 
ular model can be distinguished from that of other models by comparing, for instance, the 
forward-backward asymmetry for the process p + p — )■ /i"^ -|- /i~ -|- X as a function of the 



dilepton invariant mass, or the muon transverse momentum distribution at the LHC 28| . 
and the same asymmetry for the process + — )■ / -|- / (/ = g, /) at International Lin- 



ear Collider j29j. At first glance, another interesting feature is that the model seems to 
indicate that the LFV and DM are closely related. It implies that when the parameters 
are appropriate to satisfy the DM requirements, the LFV is relatively close to the present 
experimental bounds. In this way, the model can be confronted by the next generation of 
LFV experiments. 
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Appendix A: The constraint equations 

Here we show the constraint equations for the scalar potential given in Eq. ([2]) plus 
the terms after the (px introduction and without the Z2 symmetry. These equations are 
obtained by considering, after the spontaneous symmetry breaking, the linear terms {t^pip) 
in the scalar potential, and the solutions to the equations t,^ = are the critical points of 
the scalar potential. 

tH = Vh {2\hV^ + AhiVI + Ah2VI + AhsiVI + Ahs2VI + AhssVI + A 

= V^, [AhiV^ + 2\,,Vl + (A12 + K2)Vl + KrVl + A'^^Vj, + K,Vl + A'^^Vj^ 
-2//?i) - V2khixVhV^^ + V^^V^^{(3viV^, + (32zV^2), 

= {Ah2V^ + (A12 + K2)Vl + 2\2,Vl + !^2,Vl + A'^^VJ, + + ^axV?. 

-2111^) - kh2xVhVI^ + K,,y^3(/3i3V^^i + /323V0J, 

t^, = V^, [AusiV^ + A\,Vl + A',,Vl + 2K,Vl + Ai^l^l^ + ^nVl + A^^^^ - 
t^, = V^, [Ahs2V^ + A',,Vl + A'22F4 + Ai2V;2 + 2K2VI + A23l^3 + AsxVj, - 2/x^2) 
= V^, [AhszV^h + A'laV-l, + + Aisl^^^^ ^ ^^^^^^ ^ 2A,3^^'3 + ^3x1?, 

+3^3XV^03V^0X - 2/i^3) + Vi,KE.,(/3i3l^^, + P2ZV^,) + \/^,\/^,(2/^i23^^3 + /Sx^^;,), 

t^^ = V^, {AhsxV^ + A\xVl + A'^xVl + Aixl^^', + A^xV^ + A^xV^ + 2A.xC 
-2tiH2xVHV^, - 2filx) - V2khixVhV^, + f3xV^,V^,V^, + (5^xVl. 
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